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Abstract: The one-loop effective potential for non-relativistic bosons with a delta function 
repulsive potential is calculated for a given chemical potential using functional methods. After 
renormalization and at zero temperature it reproduces the standard ground state energy and 
pressure as function of the particle density. At finite temperatures it is found necessary to include 
ring corrections to the one-loop result in order to satisfy the Goldstone theorem. It is natural 
to introduce an effective chemical potential directly related to the order parameter and which 
uniformly decreases with increasing temperatures. This is in contrast to the the ordinary chemical 
potential which peaks at the critical temperature. The resulting thermodynamics in the condensed 
phase at very low temperatures is found to be the same as in the Bogoliubov approximation 
where the degrees of freedom are given by the Goldstone bosons. At higher temperatures the ring 
corrections dominate and result in a critical temperature unaffected by the interaction. 



1 Introduction 



Bose-Einstein condensation is central to much of our understanding of phenomena 
in condensed matter physics It is one of the simplest processes where quan- 
tum effects manifest themselves on the macroscopic level when a finite fraction of 
the non-interacting bosons in a system start to occupy the lowest energy level. Al- 
though all the particles will then be in the same quantum state at zero temperature, 
this condensate does not have real long-range order and does not truly represent 
a different phase. It was Bogoliubov |2j who first showed that a short-range re- 
pulsion between the particles is necessary in order to have a real condensate with 
the particles in a new physical phase which is superfluid. His description of the 
condensation of interacting bosons has since then formed the basis for a much more 
detailed understanding of these important phenomena ||. 

Until very recently the only physical Bose-Einstein system of non-relativistic 
particles exhibiting a phase transition at low temperatures, was liquid He 4 . But 
here the particle density is so high that it is a strongly interacting system so that 
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perturbation theory around the free system does not work |f|. However, with the 
recent experimental progress made in connection with magnetically trapped bosons 
in the gas phase [[|, the situation has radically changed and systems of weakly 
interacting bosons can now be studied. These were theoretically investigated in a 
series of papers by Lee, Yang and their collaborators forty years ago using methods 
from statistical mechanics [5], |], [7], ||. Their many results still represent to a large 
degree the state of theoretical understanding of weakly interacting boson gases. 

In the normal phase one can assign a definite number of particles to the system 
while this is impossible after the condensate has formed 0. Bose-Einstein conden- 
sation of interacting particles is therefore the oldest and still probably the simplest 
example of spontaneous symmetry breakdown which today lies at the very heart 
of modern elementary particle theory |I0|| . Since these theories are relativistic, it 
is not obvious how Bogoliubov's method can be used in this case. Instead one has 
developed very powerful methods based on Feynman's path integral formulation of 
quantum field theories [II] which allow the calculation of the corresponding effective 



potentials in a very systematic way [12, 13, 14 



This approach to spontaneous symmetry breakdown based upon functional me- 
thods has not yet been used to the same extent in the study of Bose-Einstein con- 
densation of non-relativistic systems although the basic elements are already in a 



modern textbook [15]. In two dimensions it has been used by Lozano||T6[. It is 
not clear if this approach is equivalent to the Bogoliubov method or not. This is 
what we have set out to investigate here. A first step in this direction was taken 
several years ago by Kapusta |17| who considered interacting systems of relativistic 
bosons at non-zero chemical potential and their condensation at low temperatures. 



More recently, Bernstein and Dodelson [IS] and Benson, Bernstein and Dodelson 



19(1 extended these relativistic calculations and also considered the non-relativistic 
limit. We will here show that their finite-temperature results are incomplete in that 
they have not included the contributions from the ring or daisy diagrams which are 



known to be essential at non-zero temperatures [0, [H], [21], ^TJ. Since then, the 



non-relativistic Bose gas has also been studied by Stoof and Bijlsma [22 



The use of functional methods in quantum statistical physics of non-relativistic 
systems has not yet as widespread as for relativistic systems. One of the best 
introductions have been given by Popov [|23|| . Our approach is different and more 
along the lines used in relativistic quantum field theories |23| , but we will to a large 



degree reproduce his results. The necessary formalism is established in the next 
section where we will derive the thermodynamics of a gas of free bosons in this 
language. We will work with the two real components of the field instead of the 
complex field itself and its conjugate which is usually done in condensed matter 
physics [25], [2~6| . We find this choice of variables especially advantageous in the case 
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of interacting particles at zero temperature considered in Section 3. We calculate 
the effective potential and free energy at non-zero chemical potential in the one-loop 
approximation where we include the quantum effects of the fluctuations around the 
classical solution. After removing the divergences in the theory by renormalization 
of the coupling constant and the chemical potential, we find the ground state energy 
of the hard-core bosons to be in exact agreement with the standard results of Lee 
and Yang |J. In Section 4 we extend the calculation of the effective potential to 
finite temperatures using the imaginary-time formalism. We discover that the one- 
loop approximation is no longer consistent with the Goldstone theorem requiring the 
excitation spectrum to be linear in the long wavelength limit. The problem is solved 
by including the so-called daisy or ring corrections to the boson propagator. These 
important contributions to the free energy of relativistic quantum field theories at 
finite temperature were first discussed by Dolan and Jackiw |L3] and Weinberg |L4 



They have recently become of importance in the connection with the standard model 
of elementary particles at finite temperatures [ED, I2TJ. In the last section we discuss 



the obtained results and compare with what has been obtained by other methods. 
This functional approach allows now in principle a systematic calculation of higher 
loop corrections to the thermodynamics of the interacting boson gas. 



2 Functional methods for the non-relativistic bo- 
son gas 

The wavefunction ip = ip(x., t) for a free, non-relativistic particle of mass m satisfies 
the Schrodinger equation 

idt^ = -T^V 2 ^ (1) 
2m 

when we use units so that h = 1. In the second quantized description of a system 
of many such particles ifj{x,t) becomes the corresponding quantum field. The wave 
equation (|1]) is then the classical equation of motion. It follows from the Schrodinger 
Lagrangian 

£ = J_| V ^| 2 (2) 

2m 

Constructing now the Hamiltonian H and the number operator N = fd 3 xip*tp, the 
grand canonical partition function for the is then given by 

E(P,ii) = Tte-^ H -'^ (3) 
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where (3 = 1/T when the system is in thermal equilibrium at temperature T and 
chemical potential fi. The Boltzmann constant is taken to be ks = 1. Rewriting 
now the trace as a path integral |pLT | a one obtains the functional integral 



H(/3, pi) = J T>ipViP*e-fo dT f d3 * c ^*^ (4) 

where the field ip = ?/>(x, r) is now function of imaginary time r = it. Its dynamics 
is governed by the Euclidean Lagrangian density 

C E = rd T ip + I VVf - (5) 
2m 

where we have included the contribution from the chemical potential. 



2.1 Complex field formalism 

Since the integral (|J) is Gaussian, it can be immediately evaluated. The result for 
the free energy Q = — In H//3 can then be written as 



tt(f3, n) = -Tr ln(«9 r - V 2 /2m - fj.) 



(6) 



The functional trace is given by eigenvalues of the indicated operator. They are 
found by expanding the Bose field in plane waves 



/ 1 00 

V>(x,r) = J— £ E^n k e ik ' x+ ^ 

\ P V n=-oo k 

where uo n = 27m/ (5 are the corresponding Matsubara frequencies. Then we have 



(7) 



^(AaO = ^E E \n{-iu n + £ k - fj) 







k n=—oo 



where £k — k 2 /2m is the single-particle energy. One can regularize the divergent 
sum by taking the derivative with respect to ek = £k — A 4 - Using then the standard 
sum 



E 



U! 2 + U! 2 OJ 
n=—oo n ' 



1 1 

2 + - 1 



and integrating back, we have 

W,^) = E[l e k + Tln(l-e- /3e k 



(9) 



(10) 
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after discarding an infinite constant. The pressure in the gas is now simply P = 
—Q/V. Taking the infinite- volume limit, it becomes 



I 



d 3 k r 



(e k - / u)+Tln(l-e- /3(e k-M) 



) 



(11) 



(2tt) 3 L2 



The first term is the zero-point energy which is the only contribution at zero temper- 
ature and is usually without physical content. However, in the next section where 
we discuss the interacting theory, it will be important. 

In order to find the equation of state, we want the pressure as a function of the 
density p = dP/dfi. Ignoring the zero-point energy in (JTTJ) , we get the ordinary 
result 



from which one can calculate the chemical potential fi as function of temperature 
and density. From the pressure one then has the equation of state. It involves a 
critical line, determined by // = 0, i.e. 



For densities above the critical value, the pressure in the gas is seen to be independent 
of the density. This condensed phase has thus an infinite compressibility which shows 
that it is non-physical. A short-ranged, repulsive potential between the particles 
will solve this problem and the Bose-Einstein condensate will become a physical 
superfluid [|TJ. 

2.2 Real field formalism 

When we later consider the interacting gas, we will find it more convenient to take 
the two real components of the complex field ip as independent field variables, 




(12) 




(13) 




(14) 



The Euclidean Lagrangian (|5D can then be written as 



(15) 
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where equal indices are summed from 1 to 2. We have here introduced the matrix 
operator 



M, 



a b 



2m 



ab 



(16) 



where e a b is the antisymmetric tensor in two dimensions with en — 1. After express- 
ing the corresponding action in terms of the complex Fourier components in (0), the 
partition function (f|) is now given by the functional integral 



E(j3,fi) = JvikVifo exp 



2 XX e lkV>]Vl + e 2 k^2^2) + U n{lPH>2 ~ ^l) 
n,k n,k 



(17) 



The two energies e lk and e 2 k are here actually both equal to e k = £ k — ji. It should 
also be clear that in this integral we have suppressed the Fourier indices on the field 
components. 

From the form of the partition function ([17]) we see that in terms of the real 
field components, the free, non-relativistic field theory is at the fundamental level 
interacting. The diagonal part of the action moves the fields ipi and ip2 at fixed time 
with the propagators 



■ D (0) 



1 
eik 



■ D (0) 



1 

e 2 k 



while the kinetic term provides the interaction with the simple vertex 
O- : Ur, O 



(19) 



of strength given by the Matsubara frequency. The full, free energy is now obtained 
in standard perturbation theory which is almost trivially solved to all orders in the 
interaction. First, we need the partition function of the free theory 



D^){D^)2 exp 



1 



XX e ik^i + e 2 k^2^2) 



= det 2(e lk )det 2(e 2k ) 
Taking the logarithm, we find the non-interacting result 



(20) 



pn 



[Tr In e lk + Tr In e 2k ] 



In 



(21) 
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where the closed loops denotes the trace over the variables in the propagators. The 
kinetic interaction will now perturb the fields in these two loop diagrams. A 1-field 
will be converted two a 2-field and vice versa with the coupling constant u n . Since 
the free propagator (V'iV^q — 0, only loops with an even number of interactions 
will contribute, i.e. with the same number of 1- and 2-fields. We then find for the 
full free energy 



pn 



iln 
2 



iln 
2 



1 

i + 2 



E 



ln(e lk e 2k ) + 



l(eike 2 k) 2(e lk e 2 k) 2 3(e lk e 2k ) 



+ 



+ ■ ■ 



EM^n + ei k e 2k ) 



'.22) 



Again regularizing as in and summing over the Matsubara frequencies with the 
help of (|9|), we recover the the standard Bose-Einstein free energy (JlQ|). 



2.3 Free propagators of real fields 

The propagators (|18D move only the fields at fixed time. Motion in time is induced 
by the kinetic interaction in (|T5|). Its full effect can easily be calculated in pertur- 
bation theory. For the 1-field, when we again consider a Fourier mode with a given 
momentum k and Matsubara energy u n , we find: 



Du = (M)= 

= + o--— o + o— — o o~— o + 

= D®+D®{u, n )D<g(-u, n )D$ 
+ D$(u n )D$(-u n )D$(u n )D$(-u n )D<$ + ■■■ 
= _ e 2k 

1 + ^ + ei k e 2k 

Similarly, we find 



(23) 



+ — o a—- + — o o~ — o a---- h 

(24) 

^ + ei k e 2k 
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and 

D 12 = (ViV>2*>= — »" 

= a---- + o--— o a---- H 

= (25) 

+ ei k e 2k 

while D21 = —D12. These results can simply be summed up in the Dyson- Schwinger 
equations for the full propagators 

D ab = D^ + D^jD db (26) 

With the free propagator ([TJD, = and the non-diagonal self energy IT^ = 
e c dU) n) the equations are easily solved to give the same results as above. 

We will later use the Schwinger- Dyson equations to calculate the propagators 
when the bosons have a short-range interaction. Here the fields are essentially free 
and the propagators can be obtained directly from the matrix operator flip]). Its 
Fourier transform is 

M ab =(j 7") (27) 
\ ^ n e 2 k J 

Taking the inverse, we then simply have 

D*=(*jk) = M£ = - ir ± ( J M (28) 

Un + eike 2 k \ ~ Un eik / 

which is seen to agree with the previous results. 



3 Hard-core bosons in the one-loop approxima- 
tion 



We will here consider the idealized case of bosons having only a repulsive interaction 
potential V(r) at short distances. The thermodynamics of the gas will then be 
mostly independent of the detailed shape of the potential which will only enter the 
results via the S-w&ve scattering length 



^jdhViv) (29) 



which is positive. This is equivalent to saying that the potential is a 5-function, i.e. 
V(r) = 2A<5(r) with the coupling constant A = 2na/m. In the second quantized 
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theory it will correspond to an interaction term A (if)* if)) 2 in the Lagrangian. While 
the coupling constant A would be dimensionless in the corresponding relativistic 
theory, it is not in the non-relativistic description we use here. Let us comment 
briefly upon this point. 

The Lagrangian of a real relativistic scalar field \l/ can be written as 

C = -d^W - - (30) 

in real time formalism. Here \x = . . . 3 and Ao is the relativistic coupling constant. 
Since the action must be dimensionless, the field takes on the same dimension as 
the inverse time or inverse distance in units where the velocity of light c — 1. The 
coupling Ao is therefore dimensionless. We now take the non-relativistic limit by 
letting m — > oo. Before doing so, we introduce the non-relativistic field if) through 

* = (e- imt i) + e imt i)*) (31) 

This leads to a number of terms in the Lagrangian which oscillate with frequency 
2m. They may be dropped as m — > oo. The resulting non-relativistic Lagrangian 
takes the form 

C = iif;*d t if; - -^|W| 2 - A(^) 2 (32) 
2m 

with a non-relativistic coupling constant A = 3A /2m 2 . This coupling is obviously 
not dimensionless. 



3.1 The classical ground state 

Including the above interaction, the Euclidean Lagrangian ([|) describing the bosons 
is changed into 

C E = rdr?P + ^-\Vip\ 2 - tiif)*if> + A(^» 2 (33) 

In the classical limit at zero temperature the system will be in the lowest energy 
state. The field will then attain a constant value given by the minimum of the 
classical potential 

U(iP) = -iiif)*if) + \(ip*ip) 2 (34) 

It is invariant under the U(l) phase transformation if)(x) — > e ie if)(x) and thus de- 
pends only on the modulus \ip\ of the field. In Fig.l the classical potential is plotted 



9 



for the two cases /i > and p < 0. We see that in the first case, we will have a 
ground state with spontaneous breakdown of the U(l) symmetry in which the field 
takes the classical value \ip\ = J p/2\. We will in later sections see that \i will be 
negative at high temperatures and the gas will be in the normal state with |?/>| = 0. 



U(y) 



ja < / 









Figure 1: The classical potential U(tp) plotted as function of the field modulus \ip\ in arbitrary 
units. 

To be more quantitative, the symmetry-broken ground state has a classical pres- 
sure given by the minimum of the potential (|34]), i.e. P{p) = /i 2 /4A. The corre- 
sponding number density is then p = p/2X and therefore P = Xp 2 . This is also 
the classical ground state energy density £{p) as follows from the thermodynamic 
Legendre transform 

£(p) = W-P(p) (35) 

In this way we have already cured one of the problems of the ideal Bose-Einstein gas, 
namely the infinite compressibility of the condensed phase. What is needed next 
is the inclusion of quantum and thermal fluctuations around this classical ground 
state. 
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3.2 Gaussian fluctuations 



Denoting the fluctuating field by x, we can write the full Bose field as 

1>{x) = y/±[v + x{x)] (36) 

where v is the constant field of the condensate. Inserting this into the Lagrangian 
(|33|), we can rewrite it as 

£e = ~\v\ 2 + ^\v\ A + \xaM ahXh + \{v -X)(X-X) + \{x ■ X? (37) 

when expressed in terms of the two real components of the complex fields x an d v. 
Terms linear in x have been dropped since they will not contribute after integration 
over space. The matrix operator M ab is now 

- /V 2 \ 

M ab = ie ab d T - I — - X(v ■ v) + n 1 5 ab + 2A v a v b (38) 

Due to the £7(1) symmetry of the system, we can choose the classical field v to be 
real. Taking the Fourier transform of the operator, it becomes 

M ab =( £ *- fi + 3Xv2 ~ U 1, 2 ) (39) 

where again e-^ = k 2 /2m. 

The grand canonical partition function @ can be now written as a functional 
integral over the real fields xi an d X2 as 

E(j3,fi) = e -m-W+&) Jv X iV X 2e-fo d TJ d3 *lh^ ab x b +c mt ] (4Q) 
where V is the volume of the system and 

Cint = Xvxiix ■ X) + j(x ■ X) 2 (41) 

gives the interactions of the fluctuating field. The contribution from the first part 
of the Lagrangian being quadratic in the fields, can be evaluated as for the free 
theory in Section 2 and will be the one-loop result. Higher loop corrections due 
to the cubic and quartic terms in C int can then be systematically calculated in 
perturbation theory. These will be especially important at finite temperatures and 
will be considered in the next section. 
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In the one-loop approximation, we keep only the quadratic part of the La- 
grangian. The integral ( f4"0| ) is then Gaussian and gives the corresponding free energy 



A n(/i, v) = -V + ±v* + In detM (42) 

From ([39|) we find the determinant detM = n^kC^n + w k) where the excitation 
energy is now 



^ = \/(e k -i2 + \v 2 )(e k -i2 + 3\v 2 ) (43) 

Summing over the Matsubara frequencies using ([]) again, we obtain for the free 
energy (H) 

\ 1 r1 . n 

(44) 



i fifo «) = + ^ 4 + i £ [~^k + T In (l - e^k 



y ^' ' 2 4 v k 

As a function of the condensate v, this is the effective potential U e ff for the system 
in the one-loop approximation. The thermodynamical free energy is the value of the 
function in its minimum and equals the negative of the pressure T). Due to the 
quantum fluctuations, this will be slightly shifted away from the classical minimum 
at v 2 = fi/X. For this value of the condensate, the dispersion relation fl43|) simplifies 
to the Bogoliubov result [Q] 

k 



cu = ^e(e + 2fi) = —^k 2 + Amfi (45) 

Here, and in the following we drop the wave number index on u. In the long- 
wavelength limit k — > it becomes linear and represents the phonon excitations. 
This is a consequence of the Goldstone theorem which requires excitations with a 
linear dispersion relation whenever a continuous symmetry is spontaneously broken 
as here 0. We will find in Section 4 that the quantum corrections will not change 
this dispersion relation except that the classical chemical potential \i is replaced by 
an effective potential p, which results from the quantum self energies in the field 
propagators. 



3.3 Renormalization 

The infinite sum (fyj) is seen to be strongly divergent in the ultraviolet limit k —>■ oo. 
A similar divergence was also found in the free theory in Section 2. We can thus 
remove the strongest divergence by subtracting the first term in (|T0|) so that we 
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recover the pressure in the non-interacting gas. Taking the infinite-volume limit we 
then have for the effective potential 

A „ r d 3 k 



i m\ ft i ^ 4 f d 3 k 1 . . / 

U e f f (v,fi,T) = --v +-v +J — —[-(uj-s + {i) + T\n(l-e 



-I3u) 



(46) 



However, the integral representing the Gaussian fluctuations at zero temperature 
is still not finite. The divergences can be removed by renormalizing the coupling 
constant A and the chemical potential fi. For this purpose we introduce the counter- 
terms 

- \(5^-3v 2 5X) X 2 1 -\(5fi~v 2 6X)xl (47) 



resulting directly from the classical potential fl34|) . Again we have ignored terms 
linear in the fluctuating field \. At zero temperature we then get for the renormalized 
effective potential 

(7 e// («,A,r = o) = -f« 2 + jt- 4 

This can now be made finite by adjusting the quantities 5 fx and 5X. The temperature- 
dependent contribution in fl46|) is finite by itself and thus the full free energy will be 
finite. 

Explicit expressions for the counter-terms can be obtained by isolating the two 
leading divergences in the integral ( fl~8"D by expanding the expression (f|3|) for u for 
large values of the momentum k. Cutting off the integral at k = A, we can then 
rewrite U e ff as 

U eff (v, fi,T = 0) = - V + - v 2 + ±5X 

1 r A d 3 k ( - AV\ 1 f A d 3 k / 2 AV\ 
+ -/ -, — r^[u-e + u-2Xv 2 H +-/ t — [2Xv 2 

2 J (2tt) 3 V p 2e J 2 J (2vr) 3 \ 2e J 

The divergences in the last integral are removed by taking the counter-terms to be 

f A d 3 k . A 



and 



., /-a d 3 k X 2 X 2 

(27T)' 3 £ 7l z 
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These give the same renormalized parameters as obtained by Benson, Bernstein 



and Dodelson |Tj|. Jackiw [^7] has also obtained the same renormalized coupling 
constant from considerations of bound states of non-relativistic particles in a 5- 
function potential. The unrenormalized coupling constant Ao is thus given in terms 
of the cut-off by 

1 = 1- A («> 

A A n z 

For fixed, renormalized coupling A it increases with the cut-off. In the corresponding 
relativistic theory, this increases is logarithmic instead of linear. 

Having removed the divergences, we can now let the cut-off A — > oo. We then 
have the finite result 

U e/f (v,p,T = 0) = -7p + ^ 4 

1 f d 3 k ( nx 2 AV\ 

oo- E + /I- 2Xv 2 H (52) 



2 J (2ir) 3 \ ^ 2e 

With definite values for the couplings fi and A it can be calculated by a numerical 
integration. As a function of the condensate value v it has the same shape as in 
Fig. 1. It will have an imaginary part in the region below the classical minimum at 
v 2 = fi/X whose physical significance has been discussed by Weinberg and Wu 



3.4 Condensate and pressure at zero temperature 

Because of the quantum fluctuations, the minimum of the effective potential (|52|) is 
shifted away from the classical value. In the new minimum the derivative 
(dU eff /dv)^ = 0. With 

=(25-2/i + 3A^ (53) 

we find the minimum to be at the condensate value 

2 ii f d 3 k f2e-2jj + 3Xv 2 , Xv 2 



The integral represents here the effects of the quantum fluctuations. To this order 
in perturbation theory we can evaluate it using the classical value v 2 = \ijX for the 
condensate. It then gives 

2 // f d 3 k (2e + /i ( /i 
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where now the Bogoliubov frequency u is given by (ffij). We will in the next section 
show that the quantum shift of the classical minimum is due to the self energies 
of the interacting fields. These will effectively change the chemical potential from 
the classical value p to a new value p which is just A times the right-hand side of 
equation (|55|) at zero temperature. 



With the above value for the condensate in the minimum of the effective poten- 
tial, we find by insertion into (|52|) the value for the thermodynamic pressure at zero 
temperature. To lowest order in the quantum correction, it can be written as 

™ - £-5/SH (56) 

The full density of particles is now given by the derivative p = dP/dp, i.e. 

'-a-s/IMH- 1 ) (57) 

Most of the particles are in the condensate with k = 0. Their density p c = v 2 /2 
follows directly from (|55|) . The density p e = p — p c of particles in excited states with 
k 7^ is thus 

lfd*kfe + p_^ 



re 2 J (2tt) 3 V to 

It is caused by the hard-core repulsion between the particles and was also obtained 
by Benson, Bernstein and Dodelson ]19 . 



With the Bogoliubov dispersion relation fl£5|) for the excitation energy u, we can 
now evaluate the pressure in (|56|). Besides elementary integrations, it involves the 
integral 

f de eyfe + 2p =\{e + 2pfl 2 -\p{e + 2pf' 2 (59) 

J Do 

By construction, we now get only a non-zero contribution from the lower limit e = 
of the integrals which gives 

„, N u 2 8m 3 / 2 ,- /9 . , 

PM = Tx ~ W ' < 6 °) 

The last term represents the one-loop, quantum corrections to the classical result 
in the first term. This agrees with the original results of Lee and Yang who 
considered the same system of interacting bosons with a hard-core repulsion within 
the framework of quantum statistical mechanics. 
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Equation (|57]) enables us to relate the chemical potential to the particle density. 
It can be obtained more directly by just taking the derivative of (BO) with respect 

to /!, 



P 



3/2 



(61) 



2A 3tt 2 

The sign of the second term here is apparently opposite to what was obtained by 



Bernstein and Dodelson [[18] . It must be negative as here in order to have an increase 
in the zero-temperature pressure because of the repulsion between the particles and 
not. By inversion, we then obtain for the chemical potential to lowest order in 
perturbation theory, 



8A 

p = 2Xp + — — (2m\p) 

6TX 



3/2 



(62) 



Similarly, from ( |58"D we find the density p e = (mp) 3//2 /37r 2 of particles in excited 
states. Expressed instead in terms of the full density, it is 



Pe 



37T 



(63) 



when the coupling constant A is replaced by the scattering length a. This is just the 
textbook result iUI. 



With the chemical potential from (|6^) 

P{p) = Ap 2 + 



we now get for the pressure (E0|) 

4 m 3/2 



5tt 2 



-(2Ap) 



5/2 



(64) 



The corresponding energy density follows then from the Legendre transformation 
£ = up — P and is 



8{p) = Ap 2 + 



2 " m3/ %oX^5/2 



157T' 



-(2XpY 



(65) 



It was first calculated by Lee and Yang || using the binary collision method and 
also by Lee, Huang and Yang who used instead the pseudopotential method. 
Usually, it is expressed in terms of the scattering length and takes then the form 



£ 



2na 



m 



15 V 7T 



(66) 



As a consistency check, we see that it reproduces the pressure (pjj) using the standard 
definition P = p 2 d(£ / p)/dp. 
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4 Effective potential at finite temperature 

Ignoring the counter-terms, the effective potential in the one-loop approximation 
and at finite temperature was obtained in ( fl6|) as function of the condensate v. The 
excitation energy uj is now given by the general formula (f43|). Taking the derivative 
with respect to v, we then find that the minimum of U e ff is shifted to 



due to thermal fluctuations. Here, is a function of temperature for fixed total par- 
ticle density p. The above equation then gives the condensate as an implicit function 
of temperature. The contribution from thermal fluctuations is finite. Therefore, the 
divergences are removed by the counter-terms already introduced at zero tempera- 
ture. 

With thermal fluctuations present, v 2 may differ considerably from /i/A. From 
(|43|) we see that the excitation energy no longer satisfies uj oc k in the long- 
wavelength limit. In other words, at finite temperature when the classical field takes 
on a modified value, the Goldstone theorem seems to be violated in the one-loop 
approximation considered up to now. 

It is easy to see how the effective potential can be improved so that the Goldstone 
theorem is restored. At finite temperature there must be additional effects taken 
into consideration which changes the thermodynamic chemical potential /i into an 
effective chemical potential fl so that the condensate is again given by v 2 = /i/A. 
This will give a linear dispersion relation at low energy. The minimum of the effective 
potential will then be at 



where Vq is the zero-temperature result (|5^). This is now consistent with the Bo- 
goliubov dispersion relation 



which will be verified in the following where the ring-improved effective potential is 
derived. 

The value of the condensate decreases with increasing temperatures and becomes 
zero when /2 = 0. This defines the critical temperature for the system. At higher 
temperatures it is in the phase of unbroken symmetry where the dispersion relation 
simplifies to u = e — p, as it is for a system of free particles. 




(67) 






(69) 
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4.1 Ring corrections to the effective potential 



The chemical potential /i represents the energy of a adding or removing a single 
particle from the system of interacting bosons. It represents the self energy of the 
complex field if). We see from (^) that in the symmetry-broken phase each of the 
two real modes of the field has self energies equal to respectively /i — 3\v 2 and 
fi — Xv 2 within the Gaussian approximation. Including the interactions in the full 
Lagrangian ( pJTD these values will be modified by radiative loop corrections. When 
the system is at a non-zero temperature, the additional contributions to the self 
energies will depend upon temperature. This will then enable us to define a new, 
temperature-dependent chemical potential Jl. 

In the Gaussian approximation used in the previous section we found that the 
one-loop contribution to the effective potential (42) followed directly from the inverse 



of the field propagator = ( XaXb ) • After Fourier transformation we know from 
(|39|) that it has the matrix form 

^ = ( e z ) <7o) 

where = — \i + 3At> 2 and e 2 k = £k — A 4 + ''■ This lowest correction to the 
classical result is thus 

Qi = ^Tr In IT 1 = ±^ln(^ + e lk e 2k ) (71) 



which gave the standard expression (flip . 

The matrix (]70 ) is the inverse of the free propagator which is used in the sim- 



plest version of the one-loop approximation. It was pointed out a long time ago in 



connection with the effective potential for scalar, relativistic theories 0, [TJ| that 
self energy corrections to the field propagator gave important contributions to the 
free energy at finite temperatures. More recently these daisy corrections have been 



investigated in more detail |2(J and are now often sometimes called ring corrections 
to the effective potential. We will here see that they also play a crucial role in the 
non-relativistic theory in saving the Goldstone theorem at non-zero temperatures. 
A technique which incorporates the use of resumed propagators in a self-consistent 
way is the effective action for composite operators |29]. In the Appendix we show 



that this method leads to the same result as found here, to the order considered. 

Denoting the interacting propagator by D a b, it will in general satisfy the Dyson- 
Schwinger equation (p6|). If n a ^ is the full, one-particle irreducible self energy, the 
equation then becomes 

D ab = D ah + D ac U cd D db (72) 
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Figure 2: A graphical representation of the Schwinger-Dyson equation (|7^). The thin and thick 
lines represent the free and interacting propagators, respectively. 



as shown in Fig. 2. Using now this propagator in the one-loop approximation, we 
get the modified contribution 



fii = — Trln/r 1 

2{3 



(73) 



to the effective potential. Since we can write (721) in the form 



d* 1 = d-, 1 - n a 



(74) 



we have 



Hi = — TrlnD-\l -DU) 



f2i - — Tr 

1 2(3 



DU+-(DUf + -(DUf + 



(75) 



These additional contributions to the lowest order result are called ring corrections 
as seen in Fig. 3. 




+ 




+ 




+ 



Figure 3: A sketch of the lowest order ring diagrams. 



Perturbative calculations will in the following show that the off-diagonal self 
energies II12 = — II21 vanish for zero external energy. This follows in general from 
time-reversal invariance, which implies that the off-diagonal self energies must be 
odd in the Matsubara frequency u n . The structure of the interacting propagator 
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will thus be the same as for the free propagator (|70|) but with the diagonal elements 
changed to 

eik = - A* + 3Af 2 - rin (76) 

and 

e 2 k = - A* + Xv 2 - n 22 (77) 

The ring-corrected one-loop contribution to the effective effective potential ( [7^ ) is 
thus 

^i = ^E ln (^ + ei k e 2k ) (78) 

^ n,k 

and will be of exactly the same form as the lowest order result ( pH| ) except for the 
modified dispersion relation 



^k = V( £ k - /i + 3A^ 2 - n u ) (e k - fx + Xv 2 - II 22 ) (79) 

for the elementary excitations above the ground state. 

The Goldstone theorem can now be satisfied also at non-zero temperatures if 
the chemical potential is related to the value of the field v in the minimum of the 
effective potential by 

H = Xv 2 - n 22 (80) 

This is our form of the Pines-Hugenholtz relation [^] which is usually written is the 
complex field basis. A simple derivation of this theorem is given in [^2j . Introducing 
the effective chemical potential 

/u = /i + n 22 (8i) 

the expectation value of the field when the system is in thermal equilibrium is then 
given by v = (/2/A) 1 ^ 2 as it is at zero temperature. This will be demonstrated in the 
next section. 

The dispersion relation ([75| ) now becomes 



^k = VWk + 2/x + n 22 - nu) (82) 

and is by construction linear in the momentum k in the long- wavelength limit. When 
the self energies are calculated in the following to the lowest order in perturbation 
theory, the difference II 22 — Iln will be negligible. The dispersion relation at finite 



20 



temperatures is thus the same as at zero temperature (fj5|) when the chemical po- 
tential \x is replaced by fl, i.e. it has the desired form d69|). Similarly, the inverse 
full propagator ( j74[) has the form 



d; 



if the Goldstone theorem is to be satisfied. 



33) 





Figure 4: Four-point vertices. Solid lines represent the -Dn-propagator and dashed lines the 
_D22-propagator. 



4.2 One-loop contributions to the self energies 

The interacting propagators and the corresponding self energies can be obtained 
from the full partition function 5 of the system. It is given by the functional integral 
(W\)- In the exponent we see that there is the usual coupling (—A) between four 
excitations \ as shown in Fig.4 but also a new coupling of magnitude (— Xv ) between 
three excitations due to the presence of the condensate ( ip } = v / \/2. Since we have 
chosen this expectation value to be real, this latter coupling will generate only the 
two vertices shown in Fig. 5. 



Figure 5: Three-point vertices due to the non-zero condensate. 



Separating out the free partition function S , we will here calculate the part 
Hi = H/S due to the interactions in lowest order of perturbation theory, i.e. include 
terms to order A. Since the chemical potential is now defined to be fl = Xv 2 , we 
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must include all diagrams where the four-excitation coupling occurs once or the 
three-excitation coupling occurs twice. Only connected diagrams will contribute to 
In Hi and they are given with the corresponding combinatorial factors in Figure 6. 
The lines in the diagrams represent here the free propagator D ab given in (70). 



In E, = 3/4 OO + 3/4 C X . C + 1/2 OC ^ 



+ 3 



9/4 



OOO + 3- ooo ♦ - oco 



9/4 ' • « . + 3/2 



OCO + 1/4 coo 



-oco -ooo —coo 



3/4 +3/4 +1/2 «CC> +3/2 





+ 3 O +1/2 €Z> +3 4 




o 



Figure 6: Two- and three-loop diagrams contributing to ln^,/. Only two-loop diagrams are 
considered here. The lines represent the free propagator D a b. 

We can now derive the full propagator D ab = ( XaXb ) from the partition function 
H by taking the derivative with respect to the free propagator D ab as one can see 
from the functional integral fl40l). One then finds 

51nS 51nS n 5 In Si „ 
Daf, = -2 t = -2 r + 2D ca D bd 84 

The first term is just the free propagator while the last term gives the self energy. 
In terms of diagrams, the derivative is obtained directly from the bubble diagrams 
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in Fig. 6 by opening up the corresponding lines in all of the loops. For example, 
Dn is given by the diagrams shown in Fig. 7. 



Ai = x x +3 x- 




18 x 4 \ x+2x 1 f x+12x- 




-O. 



+3 x o- --»*--»- 



+4 x — o- 



+12 x- 



hQ-h— x +4x-^>-Q x + 4x Q- 



Figure 7: The lowest order diagrams contributing to the interacting propagator Dn. Crosses 
indicate that the external propagators should be included. The lines represent the free propagator 

Dab- 

By cutting off the external lines in these diagrams, we can read off the different 
self energies U c d shown in Fig. 8. This is equivalent to using ([72]) in lowest order of 
perturbation theory where it gives D ab = D ab + D ac Yl cd D dh from which Ii cd can be 
isolated. 

An infinite set of higher loop contributions to the free energy can now be ob- 
tained by replacing the free propagator D ab in the diagrams in Fig. 6 with the full 
propagator D^. This iteration then generates diagrams which correspond to the "su- 
per daisies" of Dolan and Jackiw [O. They have more recently been investigated 
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Fi gure 8: Self energies to lowest order in the loop expansion. The external propagators should 
not be included. 



in more detail by Carrington and others In this improved perturbation 

theory the self energies in Fig. 8 are then to be calculated with the full propagators. 
In order for the Goldstone theorem to be satisfied, we know that D ao must have 
the form ([83]) involving the effective chemical potential pL. If this approximation 
method is consistent, the unknown function fi(T) can then be determined from the 
Pines- Hugenholtz relation fl81"|) which now will constitute a gap equation. 



There are other diagrams in addition to those given in Fig. 6 and Fig. 7, but they 
all vanish. Tadpole diagrams, such as — Q , are zero because (xi) = = ( % 2 )■ 
Diagrams with a 5 12 -loop -Q> are odd in the Matsubara frequency and thus vanish 
upon summation. The perturbative contributions to the diagonal self energies are 
non-zero. They depend in general on external momentum and energy. We have 
assumed a spatially constant condensate and can thus take the external momenta 
in the self energies to be zero. For the external energies we will also take these 
to be zero. This is obviously an approximation. At sufficiently high temperature, 
the zero-energy Matsubara modes should dominate the partition function and the 
approximation should then be reasonable. Especially for II 22 which we will need in 
the Pines-Hugenholtz relation (^1|), the approximation should be good because the 
pole in the corresponding real-time propagator is at zero energy. The off-diagonal 
self energies n 12 and n 2 i in Fig. 8 both vanish at zero external energy. As earlier 
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stated, this follows to all orders in perturbation theory from time-reversal invariance. 



From the diagrams in Fig. 8 we then find the self energy 



n 



11 



1 f ^ 



3(-X)D 11 (k) + (-X)D 22 (k) + 18(-Xv) 2 D 11 (k)D 11 {-k) 



+ 2{-Xv) 2 D 22 {k)D 22 {-k) + \2{-Xv) 2 D l2 {k)D 12 {-k) 



i5) 



where we use k = (oo n , k). With the full propagators from (^) this can be simplified 
to 



nn 



A 



d 3 k 



4e + 2jj 2Xv 2 (10e 2 + 4ep + A/2 2 -Quj^) 



uj 2 + uI 



(86) 



(uj 2 + uj 2 n ) 2 

where uo 2 = e(e + 2/2). The Matsubara summations are now easy to do. From the 
standard sum 

t 

UJ 



E 



follows directly by differentiation 



1 1 

2 + - 1 



7) 



E 



10 



(u 2 + ulf 



2uo 



1 1 

2 + e^ 7 



- 1) 



With co>^ instead of uo 2 in the numerator on the left-hand side, we get the same 
result except for an opposite sign in front of the last term. The zero-temperature 
self energy is then found to be 



n n (r = o) 



-A 



-A 



d 3 k 
(2^)3 

d 3 k 
(2^)3 



2e + 4/2 fx(5e 2 + 2eft + 2/2 



r.2\ 



UJ UJ 

2e + p, 2p,(e — /2) 2 



39) 



In the last line we have regrouped the terms to show the difference between rin and 
U 22 found below. The leading divergences will be cancelled by the counter-terms. 
This will explicitly be demonstrated in the case U 22 which we evaluate next. From 
the diagrams in Fig. 8 it is given as 

d 3 k 



n 



22 







E 



(2tt)* 



3(-X)D 22 (k) + (-X)Dn(k) 



+ Ai-XvfD^D^-k) +4(-Xv) 2 D 12 (k)D 12 (-k) 



X 

1 
X 

1 



E 

n 

E 



tot 



d 3 k 
(2tt) 3 

d 3 k 4e + 2/2 
(2vr) 3 uj 2 + ujI 



As + 6/i A\i{e 2 + 2e[i + uo. 



{uo 2 + uolf 



(90) 
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A summation over Matsubara frequencies gives at non-zero temperature 



n 22 (T) = -\J 



d 3 k 4e + 2/2 
(2tt) 3 uj 



1 



1 

2 + - 1 



(91) 



The first part of the integral which represents the zero-temperature self energy, is 



ultraviolet divergent. However, from (f47 ) we see that we should add in the counter- 
terms 5 [A — v 2 5X. Taking these from (|4"5| ) and (0), with the replacement /t — > ft, we 
obtain the renormalized and finite result 



n 22 (T) = -a/ 



2e + /i /i 45 + 2/i 1 

uj e uj e j3u — 1 



(92) 



Similarly, an expansion for high momenta of the zero-temperature part of fin in 
(|S9|) shows that the counter-terms 6fi — 3v 2 5X introduced in (|47|), again with fj, — > /2, 
exactly cancel the ultraviolet divergences in IIn. 

The assumption that the Goldstone theorem is fulfilled, i.e. the Pines-Hugenholtz 
relation fl30|), can now be shown to equivalent to the minimalization criterion (|68|) 
for a one-loop potential calculated with the full propagator D^. Combining (|55|) 
with fi — > /2 in the one-loop term and (p3) we find 



fi r d 3 k 
A ~ J (2^)3 

^ + In 22 (T) 



2e + /i /x 

Co" £ 



4e + 2/2 1 



e /Jw _ 1 



(93) 



when making use of the result (|92|). This is the same as equation ( jSCf ) which embodies 
the Goldstone theorem . 

As mentioned at the beginning of this section, the self energies are calculated 
using the dispersion relation uj ■■ 



'e(e + 2/2), i.e. neglecting terms of 0(1111 — Il 22 ). 
Using the renormalized self energies, we see that the difference at zero temperature 
is 



n 



ii 



n 



22 



2A/t 



d 3 k 



(94) 



The integral is finite at high momenta due to the counter-term, but is infrared 
divergent. This is caused by the exchange diagram in II n with two Z) 22 -lines which 
diverges as external energy and momenta are taken to zero. The infrared divergence 
signals the onset of long-distance effects which are not properly handled by the 
present one-loop approximation. As discussed by Kapusta |24j], one can cure the 
divergence using a non-zero external energy. However, this is not a problem here. 
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At low temperature, the self energies are of (9(A/i 3//2 ) and will be neglected all 
together, being small compared to fi. On the other hand, at high temperature, i.e. 
near the critical temperature T c , H22 is comparable to \i and the self energies must 
be included. In this regime we can use the high-temperature limits of ITn and n 22 
which below are shown to be equal to leading order. Thus, in the temperature range 
where it should be considered, Iln — H22 = 0. 

At high temperatures T ^> p, the exchange diagrams may be neglected. The 
dominating terms in the remaining diagrams go like e/{u(e^ w — 1)), and the combi- 
natorial prefactors are the same for the two self energies. Expanding (|86f ) and fl92|) 
in powers of fl/T we find to leading order 

(mT\ 3//2 

IT 22 (T) = n n (T) = -4AC(3/2) (— J (95) 

This result will be found if one uses free propagators D a ^ in the self energy loops 
since p can be neglected to leading order. 

In the opposite limit T <C p, we can take p approximately equal to fi. Diagrams 
with a Z) n -loop are suppressed by a factor T 2 //z 2 compared to the others. We then 
obtain for the leading terms of the self energy 

,rr* 10A / XT 2 m 3/2 

n»m = - («) 

The first part is seen to be of the order 0(\Zpd d ) smaller than the chemical potential 
and the second term even smaller by a factor (T/jl) 2 . To leading order, we can thus 
neglect the self energies at these low temperatures and set p, = //. Thus, the effects 
of ring corrections on the thermodynamics of the system will only appear at high 
temperature. 



4.3 Condensate and pressure at finite temperature 

Including the effects of ring corrections, we can now write down the renormalized 
finite temperature effective potential 

+ T /(^ h '( 1 -^" ) < 97) 
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From (|93D we know that it has a minimum at v 
then have the pressure 



P(p,T) 



-U eff (ji,T) = —(fi 2 



ul 2 ) 



(p/X) 1 / 2 . In the minimum we 

d 3 k 



(2tt)' 
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T 



d 3 k 
Wf 



ln(l 



-j9un 



(9* 



The total density of particles is again found by using p = dPj dp which gives 



JL 
2A 



u 1 
— + - 
2A 2 



d 3 k 
(27p 

d 3 k 
(2^)3 



- + - 

a; £ 

e + p 

UJ 



2s 



u{e^ - 1 
2(e + p) 



oj(e 
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(99) 



while the density of particles in the condensate is 
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~2 



Pc 

The difference 
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UJ 
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4e + 2p 
u(eP" - 1) 



(100) 



1 



- 1 



(101) 



represents the density of particles in excited states with k / 0. Except for the 
replacement of p with p, this result agrees with what was obtained in [|ll|. Now both 
the hard-core repulsion and thermal fluctuations cause the excitation of particles 
from the condensate. With increasing temperatures more and more particles are in 
excited states and at a critical temperature where p = 0, the condensate becomes 
zero. We then have a phase transition to the normal phase. In the following we will 
see that this critical temperature is the same as for a free Bose gas to the accuracy 
we are working. 

Equation ( p9|) gives upon inversion the chemical potential as a function of the 
particle density. The first integral on the right-hand side was done in Section 3.4, 
and we can now simply replace p with p in Eq. (59). The second integral 

d 3 k e 



(2n) 3 uj(e^ 



(102) 



must be done numerically. However, in the two important temperature ranges we 
can easily find good analytical approximations. 
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At low temperature we again take fi ~ \i. The dominant contributions to the 
integral comes from small values of e and we can thus take e + 2/x ~ 2/i in the 
denominator. Then we have the result 



7Tm 3 / 2 T 



Since the loop corrections are small, we can set \x = 2X P on the right hand side of 
( j99|) . Inversion then trivially gives 

fj, = 2 \p + ^-(2mX P f' 2 + A7T ^ 3/ ' T 9 4 (104) 
P H 3vr 2V hj 30(2Ap) 5 /2 v ; 

For comparison, the expression for the effective chemical potential takes the form 

2A , x Am 3 / 2 T 2 , , 

, = 2X Pc = 2X P --(2rnX P r-^ w (105) 

Thus, at low T, the chemical potential increases with temperature, while the effec- 
tive chemical potential, and thus the condensate density decreases with temperature. 
This temperature dependence in the same as obtained in the Bogoliubov approxi- 
mation [25|. The pressure now follows from Eq. (|98|) . As already explained, we can 
neglect the contribution coming from the self-energy. Thus 

where the temperature dependence is in the function 

evaluated in the same approximation as above. The pressure at low temperatures 
is therefore 

s ^ 2 4m 3 / 2 , ,-,„ m 3 / 2 7r 2 T 4 , . 

^ = Xp + + HaT Wp^ (108) 

At non-zero temperature we can obtain the energy density from the thermodynamic 
relation 

£ = (109) 



which now gives 



2na 9 

8 = p 2 
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i + ™J*L 

15 V 7T 



m 3/2 T 4 

+ 240^7V^ (110) 
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It is seen that the effect of non-zero temperature comes in as a T 4 contribution as 
in the Stefan-Boltzmann law for photons. In that case the excitations are mass-less 
because of gauge invariance, while here they are Goldstone bosons resulting from a 
broken, continuous symmetry. 

In the opposite temperature range near the critical temperature Tc, we must 
include the effects of the self energies. They have previously been obtained in 
which gives for the effective chemical potential (|8TD 



fi = H- 4AC(3/2) (—J (111) 



This quantity also enters the total density of particles in (|99| ) where we now to leading 
order can set fi = 0. The temperature-dependent integral ( |102j ) then simplifies to 

It = C(m (-^r) (H2) 

Since the chemical potential is now fi = 2A(p + It), we have the relation 

p = 4Ap-/i (113) 

for temperatures just below the T c . The critical temperature is defined by p — 0. 
Here the chemical potential takes the value fic — 4Ap, and hence ( |111|) gives 



r **( * V' 3 **( p V /3 (114) 

which is just the textbook result [|I|]. As seen in Fig. 9, the thermodynamic chemical 
potential /x increases continuously from zero to \ic at the critical temperature. On 
the other hand, the effective chemical potential p, decreases smoothly from zero 
temperature and is by definition zero at Te- 
lia the general formula (98) we must now keep the contributions from the self 
energies when we calculate the pressure for temperatures just below the critical 
temperature. Thus 

p = ^ 2 - ny + ) 3/2 t 5 / 2 c(5/2) - (^) 3/2 C(3/2)(a* + n 22 ) + 0(m 2 IH5) 



with n 22 given in (|95|). In order to compare with results in the literature, we write 
the pressure as a function of the particle density 

P= (£) 2 C(3/2) 2 aT 3 + (^) 3/2 C(5/2)T 5 / 2 (116) 
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Figure 9: The thermodynamic and effective chemical potentials, fx and p plotted as functions 
of temperature for g = 0.1 and m — 1. The results were obtained by matching the low and high 
temperature expressions for the chemical potentials. 



Comparing with the result of Lee and Yang || , we find that the second term has the 
wrong sign. The reason is the following: An expansion of the ring-corrected potential 
on powers of A shows that the two-loop contribution is counted twice. Higher order 
diagrams are reproduced correctly. We must correct for this double counting by 
subtracting the two-loop contribution, as given in perturbation theory, from the 
obtained result. This contribution is given by the three uppermost diagrams in Fig. 
6. Near T c the propagators are equal to lowest order, Dn = D 2 2 = £/(£ 2 + 
The total contribution is — 2(m/27r) 2 C(3/2) 2 aT 3 . Subtracting this contribution, we 
thus have 



P= ~ap 2 + (^) 2 C(3/2) 2 aT 3 + (^) 3/2 C(5/2)T 5 / 2 
m \2tt / \2tt / 



(117) 



which now agrees with 0. This result also comes out in the mean field approx- 
imation fl], ^TJ. The equation for the chemical potential is not modified by this 
correction at the present order of accuracy. In the same approximation we also find 
the the energy density 



t = — ap 

m 



&(3/2) VT 3 / 2 + 2(^) 2 C(3/2) 2 aT 3 + |(^) 3 / 2 C(5/2)T 5 / 2 (118) 
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At the critical temperature these expressions reduce to 

P{T C ) = ^C(3/2) 5 / 3 C(5/2)p 2 / 3 + (H9) 
m m 

and 

8{T ) = -C(3/2) 5 / 3 C(5/2)p 2 / 3 + -ap 2 (120) 
m m 

which in the following are shown to equal the limits taken in the normal state, 
securing continuity at the critical temperature. 

The special treatment at two-loop order could of course have been introduced 
at an earlier stage, but there would be no effects on the results found so far. When 
minimalizing the effective ring-corrected potential we discarded contributions from 
dU/dv, being of 0(X 2 ). Since the contribution from differentiation of the two-loop 
potential is of the same order, this must be discarded as well, and the minimaliza- 
tion criterion remains unchanged. In the Appendix we briefly discuss how this is 
incorporated in the effective action for composite operators. 



4.4 Equation of state in the normal phase 

At temperatures above Tc, the condensate density is zero and the effective poten- 
tial has its only minimum at v = 0. The modified dispersion relation then reads 
uj = e — p, with effective chemical potential p = /i + H22 < 0. Since v = 0, the 
exchange diagrams vanish and the expressions for the self energies are 

n„=D- s n = ^/(gN- 

r d 3 k 1 
J (2tt) 3 e P" - 1 

The Bose-integral is standard and the self energy becomes 

IT=-4A(^_) Li 3/2 (e^) (122) 
We have here introduced the polylogarithmic function 

00 k 

k=l K 

with Li n (l) = C{n). Near Tc we may approximate the exponential with unity and 
at T c the result Q95D is reproduced. To the present accuracy we may set p, — fi in 



\)D 22 {k) + {-\)D 11 {k) 



+ c.t. 
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the self energy. With vanishing condensate the classical and zero point terms do not 
contribute to the pressure, which simplifies to 



P = ' T 5 / 2 Li 5/2 (z) + 2A 



3/2 



fmT 



7T 



U l/2(z) 



(124) 



with z = exp(/3//) and z = exp(/?/i). To lowest order in A the pressure can be written 



m \ 3 / 2 
2W 



\ 3/2 

] T 5 / 2 Li 



5/2 1 



2A 



fmT\' 



Li 3/ 2 (^) 



(125) 



where the last term is the two-loop contribution. We have here used the following 
property of the polylogarithmic functions: 



^ Li n(aO = ^Lin-i(x) 



(126) 



Expanding in powers of f3p, we find that this agrees with the result of Popov P5 
Again we want to write the pressure in terms of the particle density which now 
becomes 



P 



fmT\ 3/2 T . , . AX 
1 —) Li 3/2 ^)+4A 



m 
2^ 



r 2 Li 1/2 (^)[Li3/ 2 (z)-Li3/ 2 (^)] (127) 



Inverting this relation, we find to lowest order in A 

/i = Tin U~f 2 (pK%) + 4Ap 



:i28i 



where Li" 1 is the inverse of Li n and At = y2n jmT is the thermal wave length. To 
the same order the self energy is II = — 4Ap. Thus, the effective chemical potential 
becomes 



/i = /i + II = Tln \A~j 2 {ph%) 



(129) 



which equals the usual chemical potential for a free gas. Let us again refer to Fig. 9 
where the chemical potentials p, and \i are plotted as function of temperature. The 
difference p — p, — 4Ap found just below the critical temperature, remains constant 
in the normal phase where the potentials have the same temperature dependence. 



As a function of density the pressure becomes 

A T p 



P = — ap + pT 

m 



2 5/2 



+ o( P 3 



The corresponding second virial coefficient reads 

B 2 - 



A| 

"25/2 



1 



mT 



7T 



(130) 



(131) 
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From the relation £ 



-d 1 1 d(3{f3P)\i3^ we similarly have for the energy density 
A 3 T p 



„ 4vr 2 3 m 
£ = —ap + -pT 
m 2 



1 - 



2 5/2 



(132) 



again in agreement with earlier results || [l], 0. At Tq these expressions coincide 
with the limits taken from sub-critical temperatures, found in the last subsection. 
Summing up, the Eqs. ( |108| ), (|117|) and (|130|) constitute the equation of state for 
the Bose gas at all temperatures. 



5 Discussion and conclusion 

The non-interacting Bose-Einstein gas in the condensed phase is not a thermody- 
namically stable system since its compressibility is infinite JTJ. Also, the phase 
transition from the normal phase is special in that the correlation length diverges 
below the transition temperature. This unphysical behavior is directly related to 
the absence of a real spontaneous breakdown of the U(l) symmetry of the system 
with the accompanying Goldstone bosons which here are the phonon excitations. 

With the introduction of a weak repulsion between the particles, the physics 
of the system is well-defined in both phases. We have here described the system 
using modern field-theoretic methods based on functional integrals. Instead of us- 
ing the standard formalism with complex fields as used in the condensed matter 
literature P5f , we find it more convenient to use real fields. The divergences in 
the loop integrals have been regulated by introducing a physical cutoff. Since these 
are proportional to powers of the cutoff, we could instead have used dimensional 
regularization where they consequently would be absent. 

We have put special emphasis on enforcing the consequences of Goldstone's 
theorem. In the real-field formalism the important Pines-Hugenholtz theorem then 
takes a slightly different form. The thermodynamics is most directly obtained from 
the effective potential whose minimum gives the free energy of the system. At 
non-zero temperatures we find that the one-loop effective potential does not give 
a consistent description of the symmetry breakdown and that it must be improved 
by adding in ring corrections in a self-consistent way. It then becomes natural to 
introduce an effective chemical potential p which acts as an order parameter, being 
positive below the critical temperature Tq and giving a non-zero condensate, and 
passing through zero at Tq- In this way we find /2 is a more important variable in 
characterizing the phases of the system than the thermodynamic chemical potential 
/i. In particular, for temperatures below Tq where the particle number concept 
looses some of its meaning, the physical significance of /i is not clear. 
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Our results are derived to lowest order in the coupling between the particles. As 
such, most of the thermodynamic results are already in the literature dating back 
to the pioneering calculations of Lee, Yang and collaborators obtained by methods 
from statistical mechanics. What is new in addition to the systematic use of the 
effective chemical potential, is basically a more coherent derivation in a framework 
which has to a large extent been developed for corresponding relativistic systems in 
high energy physics. There is nothing preventing us in extending the calculations 
to higher orders in the interactions. One is then forced to consider the theory as an 
effective field theory which is non-renormalizable at very high energies, the coupling 
constant having dimensions of a length. As recently shown by Braaten and Nieto 



32], the zero-temperature ground state energy to next order in the coupling constant 
can then be obtained in a most direct way using the renormalization group for this 
effective field theory. 

With the recent introduction of magnetic traps, the experimental study of the 
thermodynamics of weakly interacting bosons have entered a new phase which will 
allow a much more detailed study of this important system. Here we have consid- 
ered the particles in an open volume where the interactions dominate, while it is the 
confining potential which has the most important role in the trapped systems ||. 
The number of particles in the system now is finite, but one can still use quantum 



field theory in the grand canonical ensemble ||33|| . Even with the large effort going 
on in the field at present, there is still a long way between experimental results and 
detailed verifications of the theoretical properties as derived for the open system 
here. 
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A Appendix 



Here we give a brief review of the effective action for composite operators introduced 
by Cornwall, Jackiw and Tomboulis p9|| . The method was originally presented in real 
time formalism, but we will use imaginary time, following |34]]. A compact notation 
with f x = J d d x and J x = J(x) is used. One starts by defining the generating 
functional 



Z[J,K] = exp(-W[J,K\) 



S[<j>) + 



Z J xy 



,.K 



xyry 



(133) 
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where the field(s) is coupled both to a local source J x and a bilocal source K xy . 
S[<f)] is the classical Euclidean action. The background field $ and the propagator 
G are defined by 

™ = (134) 

OJx 

y$£=fr.*. + G m ] (135) 

xy 

A double Legendre transform then produces the effective action 

r[$, G] = W[J, K] - [ $ x J x -If [<S> x K xy $ y + G xy K yx ] (136) 

Jx A J xy 



where J and K are eliminated using (|134|) and (|135|) . The sources are then related 
to r[$, G] via the equations 

In order to establish a series expansion for the effective action one must shift the field 
in the classical action S[4>] — > S [$ + </>]. From this the free propagator in background 
field 



itf (•> - ^ 



J -Hr% 10=0 



(139) 



is obtained. A careful rearrangement of terms leads to the series [29 



r[$,G] = S[<5>] + -TrlnDoG -1 + -Tr^G - 1] + T 2 [$, G] (140) 

r 2 [$, G] is the sum of all two-particle-irreducible graphs obtained from the interac- 
tion terms in S[<& + 0] with propagator G and -Do = D($ = 0) The conventional 
effective action is now obtained from the equations 

r[$] = r[$,G ] (i4i) 

5T[<f>,G] 



5G 



g=g = (142) 
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As seen from (|138[) this is equivalent to letting K = 0. 



We now want to compare the effective action ( |140| ) to the results obtained in 
the previous sections. $ is then the background field (v , 0) and = (xuX2)- We 
only need V to O(X). Using G^ 1 = D^ 1 — U we find for the various terms 

^TilD^G - 1] = ^TiDU + £>(A 2 ) (143) 



- Tr In G" 1 = -Tr[ln IT 1 -£>![ + •■■] (144) 

Hence, to this order, the effect of the third term in ( |140| ) is to cancel the O(X) 
contribution, in terms of the propagator D in the one-loop term. To the same order 
T 2 is simply the sum of two-loop graphs with propagator G. Omitting a constant 
due to Tr In D we may to this order write 

T[v, n] = S[v] + ^Tr lnCT 1 - ^Tr^n + two-loop (145) 

If we now make the ansatz for G — [D^ 1 — U]^ 1 , as we did assuming the Goldstone 
theorem, this is seen to equal our result for the ring-corrected effective potential. 
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